Impact of COVID-19 on routine malaria indicators in rural Uganda: an interrupted time series analysis

Background In March 2020, the government of Uganda implemented a strict lockdown policy in response to the COVID-19 pandemic. Interrupted time series analysis (ITSA) was performed to assess whether major changes in outpatient attendance, malaria burden, and case management occurred after the onset of the COVID-19 epidemic in rural Uganda. Methods Individual level data from all outpatient visits collected from April 2017 to March 2021 at 17 facilities were analysed. Outcomes included total outpatient visits, malaria cases, non-malarial visits, proportion of patients with suspected malaria, proportion of patients tested using rapid diagnostic tests (RDTs), and proportion of malaria cases prescribed artemether-lumefantrine (AL). Poisson regression with generalized estimating equations and fractional regression was used to model count and proportion outcomes, respectively. Pre-COVID trends (April 2017-March 2020) were used to predict the’expected’ trend in the absence of COVID-19 introduction. Effects of COVID-19 were estimated over two six-month COVID-19 time periods (April 2020-September 2020 and October 2020–March 2021) by dividing observed values by expected values, and expressed as ratios. Results A total of 1,442,737 outpatient visits were recorded. Malaria was suspected in 55.3% of visits and 98.8% of these had a malaria diagnostic test performed. ITSA showed no differences between observed and expected total outpatient visits, malaria cases, non-malarial visits, or proportion of visits with suspected malaria after COVID-19 onset. However, in the second six months of the COVID-19 time period, there was a smaller mean proportion of patients tested with RDTs compared to expected (relative prevalence ratio (RPR) = 0.87, CI (0.78–0.97)) and a smaller mean proportion of malaria cases prescribed AL (RPR = 0.94, CI (0.90–0.99)). Conclusions In the first year after the COVID-19 pandemic arrived in Uganda, there were no major effects on malaria disease burden and indicators of case management at these 17 rural health facilities, except for a modest decrease in the proportion of RDTs used for malaria diagnosis and the mean proportion of malaria cases prescribed AL in the second half of the COVID-19 pandemic year. Continued surveillance will be essential to monitor for changes in trends in malaria indicators so that Uganda can quickly and flexibly respond to challenges imposed by COVID-19. Supplementary Information The online version contains supplementary material available at 10.1186/s12936-021-04018-0.

. At the start of the COVID-19 pandemic, there was concern that progress in the fight against malaria would be reversed due to interruption of malaria control interventions and overwhelmed health care systems, with some modelling studies suggesting that malaria morbidity and /or mortality may double due to COVID-19 [2][3][4]. In Uganda, malaria is endemic in over 95% of the country and the leading cause of morbidity and mortality, accounting for 30-50% of outpatient visits, 15-20% of all hospital admissions, and up to 20% of all hospital deaths [5]. While successes have been registered in malaria control in the country in the last decade, including increased coverage of control interventions and reductions in the overall disease burden [6,7], there is concern that these achievements may be reversed by the global COVID-19 pandemic.
The World Health Organization (WHO) declared COVID-19 a global pandemic on 11 March 2020 [8]. The Government of Uganda swiftly implemented strict restrictions and diverted personnel and resources to minimize the spread of SARS-CoV-2 and mitigate the economic burden of the pandemic. On 18 March, 2020, mass gatherings were suspended, and a 14-day quarantine was imposed on all travellers arriving in Uganda [8]. When the country registered its first confirmed case on 21 st March 2020, additional restrictions were implemented, including: (1) closure of borders except for cargo and goods on 21 March, 2020; (2) suspension of public transport and restrictions on movement of private vehicle on 25 March, 2020; (3) mandatory testing of truck drivers on 10 April, 2021; (4) a national lockdown and curfew from 19.00 to 06.30 on 30 March, 2020, initially for 14 days but eventually extended to 26 May, 2020, and (5) closure of all schools [8]. Despite these measures, Uganda's COVID-19 cases progressively increased over the period with the country registering its first 100 confirmed cases on 6 May, 2020, first death on 23 July, 2020, first 1000 cases on 9 June, 2020 and first wave peak in December 2020.
Although lockdowns have been shown to minimize the spread of COVID-19 [9][10][11], they may have major repercussions on malaria treatment and prevention including decreased access to health care, interruption of service delivery, and disruption of delivery of malaria control interventions such as long-lasting insecticide treated net (LLIN) or indoor residual spraying (IRS) campaigns. In addition, an increase in the incidence of fever cases due to COVID-19 may lead to an increase in the number of suspected malaria cases presenting at health facilities, increasing demand for malaria diagnostics, and potentially leading to over diagnosis and overtreatment of malaria and/or a shortage of diagnostics. Indeed, a modelling study by WHO predicted that in the 'worst case' scenario where LLINs campaigns were cancelled and access to effective antimalarial treatment was severely disrupted, malaria deaths could increase by > 200% in Uganda as result of the COVID-19 pandemic [12].
Interrupted time series analysis (ITSA) is a quasiexperimental study design that can be used to estimate causal effects using observational data. It is often used to evaluate the effectiveness of population-level health interventions that have been implemented at a clearly defined point in time [13][14][15][16]. ITSA has also been used to evaluate the effect of an unplanned event on health outcomes, such as the effect of the late 2000s financial crisis on suicides in Spain [17]. A key assumption of ITSA is that the pre-intervention trend would provide a suitable estimate of the post-intervention trend had the intervention never occurred. In this study, interrupted time series analysis was used to assess the effect of the COVID-19 pandemic in Uganda on malaria disease burden and case management at 17 public health facilities across rural Uganda using three years of pre-COVID data to model the baseline trend. Results from the study will improve our understanding of the effect of COVID-19 on malaria care in Uganda and in other countries with similar settings.

Health-facility based surveillance
The study utilized data collected at 17 Malaria Reference Centres (MRCs) located in rural settings from April 2017 through March 2021 (Fig. 1). MRCs are high-volume, level III/IV, public health facilities located throughout Uganda in areas with varying malaria transmission intensities. MRCs were first established in 2006 by the Uganda Malaria Surveillance Project (UMSP), a project led by Makerere University, University of California, San Francisco (UCSF), and Infectious Diseases Research Collaboration (IDRC) groups, in collaboration with the Uganda Ministry of Health. Details of the MRCs have been described elsewhere [18] but briefly, there are currently 70 MRCs collecting high quality, individual level data from all patients presenting at the outpatient clinics of the facilities. Of the 70 MRCs, 17 were selected for this analysis based on the following criteria: (1) MRC was opened and data collection commenced before January 2017; and (2) MRC had less than 5% missing data on four of the variables (suspected malaria case, malaria test done, malaria test result recorded, type of malaria test done).
At each MRC, patient data are recorded in standardized health management information system (HMIS) registers (HMIS 002-outpatient register) at the outpatient clinics. Data are transcribed from the registers into Microsoft Access databases by on-site records officers. Data collected include patient demographics, village of residence, history of fever, whether a malaria diagnostic test was performed, type of malaria test done (rapid diagnostic test (RDT) vs microscopy), results of laboratory tests, diagnoses given, and treatments prescribed. Data are not collected on whether a patient receives or is adherent to the prescribed treatment. On average, over 1,400 patient visits and 1,000 suspected malaria cases tested are captured monthly at each MRC. Adherence to malaria case management guidelines is emphasized and supported at the MRCs through training, onsite mentorship and support from regional surveillance assistants, regular continuous medical education sessions, and monthly data analysis and feedback sessions for on-site health workers. Routine data quality assessments are undertaken at each site to ensure collection of high quality data in terms of completeness and accuracy.

Laboratory procedures
All patients suspected to have malaria are routinely sent to an onsite laboratory for confirmatory testing using either malaria microscopy or malaria RDT kits. Malaria microscopy at the MRCs is conducted by skilled laboratory personnel and generally reserved for when RDTs are not available or of limited stock and upon request from a clinician. Quality control for malaria microscopy is conducted monthly at each site by a laboratory team who are part of the core UMSP team. They also provide onsite mentorship for site laboratory personnel to ensure high quality malaria testing.

Malaria control interventions
Malaria control activities at the MRCs during the study period included: (1) malaria case management as per the Uganda National Treatment Guidelines [19]   (2) promotion of intermittent preventive treatment during pregnancy; and 3) universal distribution of free LLINs. Two national LLIN distribution campaigns were conducted during the study period, in 2017-18 and 2020-2021 ( Figs. 1 and 2). The status of IRS campaigns varied across the sites included in this study as shown in Fig. 2; three sites were in districts that never received IRS; IRS was discontinued in 2017 in districts where nine sites were located, and the remaining five sites were in districts that had sustained IRS campaigns throughout the study period.

Statistical analysis
ITSA, using segmented regression, were conducted to estimate the impact of the COVID-19 epidemic on key malaria indicators. Data from all 17 sites from April 2017 through the end of March 2021 were included in the analysis. The study period included three years pre-COVID and one year after the arrival of COVID-19 in Uganda, with the start of the COVID-19 time period defined as 1 April, 2020 (Fig. 2). Outcomes for the ITSA included number of total outpatient visits, number of non-malarial visits (visits in which malaria was not diagnosed), malaria visits (visits in which a laboratory-confirmed diagnosis of malaria was made, synonymous with malaria cases), proportion of all patient visits with suspected malaria, proportion of patients tested using RDTs, and proportion of patients with confirmed malaria prescribed AL.
In the segmented regression model, the baseline trend during the pre-COVID-19 period and the slope change following the arrival of COVID-19 in the country are modelled. To estimate the counterfactual (COVID-19 time period) trend, the slope change is set to 0. The generalized segmented regression model is as follows: where Y ht represents the outcome recorded at the health facility h in month t, β 0 is the intercept at the start of the study (t = 0), T is a linear term indicating the time in months since the start of the study observation period (models the pre-COVID trend), P t is a linear term indicating the months since the start of the COVID-19 time period (models the observed change in trend after COVID onset), C t is a linear term indicating calendar month fixed effects (e.g., January, February, etc.) to model seasonality), and R ht is a vector of monthly rainfall data (mm) lagged by one month and averaged across the district-level. Poisson regression was used to model count outcomes using a generalized estimating equation (to account for clustered observations between sites) with an autoregressive of order one (ar1) correlation structure (to account for serial autocorrelation of error terms of outcomes of adjacent time periods within-sites). Fractional regression was used to model proportional outcomes.
For ITSA models that included malaria outcomes, additional linear terms representing the number of months since LLIN distribution and the round of IRS were included to adjust for time-varying confounding. Based on surveillance data from these sites, LLIN campaigns were assumed to confer 24 months of protection, whereas IRS campaigns (which used Actellic or Sumishield-based insecticides) were assumed to provide up to a year of protection. For Poisson models, the underlying population was assumed constant throughout the study period and thus no offset term was included.
Models were then used to predict the unobserved counterfactual values (i.e., 'expected' values in absence of COVID-19) for each health facility in each month (i.e., monthly number of events or proportions had the epidemic not occurred) by setting P t to zero. For count outcomes, incidence rate ratios (IRR) were calculated by summing the outcomes across the first and second six months after COVID-19 started in Uganda and dividing the observed by the expected values. For proportion outcomes, relative prevalence ratios (RPR) were calculated by averaging the monthly proportions across the two six-month COVID-19 time periods and dividing the observed by expected values. 95% confidence intervals were obtained using a 1,000 block, percentile-bootstrapping procedure, where sites were resampled with replacement.  Table 1.

Trends in disease burden and case management over the study period
Over the four-year study period, a single annual peak was observed in the number of monthly outpatient visits, with the highest number of visits observed during midyear of 2019 (> 2,500 visits, Fig. 3A). The monthly trend in the number of outpatient visits for non-malarial illness was declining at a steady pace before the COVID-19 epidemic (Fig. 3B). The number of malaria visits followed a similar trend to that of the total outpatient visits, with peaks in these occurring simultaneously (Fig. 3C). Similar trends were observed with the proportion of suspected malaria cases (Fig. 3D). The proportion of patients tested with RDT each month remained consistently high (> 70%) and almost constant throughout the study period, except in 2018 when the proportion reached > 90% and in 2021 when the proportion fell below 70% (Fig. 3E). The proportion of confirmed malaria cases prescribed AL (Fig. 3F) was consistently high, remaining > 90% from late 2019 through 2020, before a decline was observed in 2021. The mean proportion of patients with suspected malaria who received a diagnostic test is not shown in Fig. 3 because it remained consistently high across all study time points (Additional file 1: Fig. S1) and therefore a formal ITSA was not performed.  AL, observed values were lower than expected in the last six months of observation based on pre-COVID trends. Effect estimates across all measured outcomes are presented in Table 2, stratified into the first and second six months of the COVID-19 time intervals. In the first six months after the onset of COVID- 19 0.90, 0.99)). This was likely due to an increase in the percentage of patients with laboratory confirmed malaria that received a prescription of dihydroartemisinin-piperaquine (DP); this percentage rose from 0.04% in the pre-COVID-19 period to 0.3% from April 2020-September 2020 and to 3.7% in the period from October 2020 through March 2021.

Discussion
This study evaluated the impact of the COVID-19 pandemic on the number of outpatient visits, malaria disease burden and case management practices between April 2020 and March 2021 using routine surveillance data collected at 17 high-volume, public, outpatient facilities in areas of varying malaria transmission across Uganda. In this study the number of outpatient visits, malaria cases, and non-malaria visits, as well as most indicators of malaria case management, did not reflect any significant changes in the first year following the onset of the COVID-19 epidemic. However, from October 2020 through March 2021, a modestly lower mean proportion of suspected malaria patients received an RDT for malaria diagnosis compared to what was expected for that period in the absence of COVID-19. Notably, this difference was not reflected by a change in the proportion of patients with suspected malaria who underwent diagnostic testing, as microscopy was available as an alternative to RDTs. A similar trend was observed in the last six months of observation for the mean proportion of malaria cases prescribed AL.
The emergence and rapid spread of COVID-19 across the world has created massive global disruptions on health systems, social services, and economic activity [20]. The disruption in health services is expected to be magnified in sub-Saharan Africa to a greater degree than in other regions due to relatively weak health service infrastructures, low clinician to population ratios, limited laboratory capacity, and a higher burden of other infectious diseases [21]. Delivery of malaria preventive measures and care services are some of the activities most likely to be affected by the epidemic given the high prevalence of malaria in the region [22]. Indeed, there are predictions that malaria cases will increase, and malariarelated mortality may nearly double, decelerating any gains attained in the last decade [2,3,12,23]. Uganda, Table 2 Estimates of the impact of COVID-19 epidemic on disease burden during the first and second six months of the COVID-19 time period a Estimates reported as the total number of visits or mean proportion during the 1-year post-policy period b For visit outcomes, this represents an incidence rate ratio; for proportions, a relative percent ratio. Assumed a constant population at-risk over the study period one of the highest malaria burden countries in sub-Saharan Africa, was predicted to be one of the countries most likely to be affected by the COVID-19 pandemic. According to a WHO modelling study, by the end of 2020 in a scenario where LLIN campaigns are not implemented and LLIN continuous distributions and access to effective anti-malarial treatment are reduced by up to 75%, malaria deaths could increase by 200% in Uganda [12]. A study by Bell et al. predicted less disastrous effects resulting from a widespread COVID-19 outbreak compared to the impact of lockdowns on malaria programmes [24]. However, contrary to those worst-case scenarios, these data suggest that despite travel restrictions and other lockdown measures, access to care and malaria case management were largely unaffected by the COVID-19 pandemic in the first year after it was confirmed in Uganda. This could be because 75% of the population in Uganda resides in rural areas where patients seeking care either walk or ride bicycles to the health facilities [25], making access to care possible even with restrictions to motor vehicle movements. In addition, the health system is structured in such a way that most of the health care staff are housed at or near the facilities; furthermore, travel permits were provided for health workers who needed to travel long distances to their duty stations, preventing the staff shortages that would likely have limited access to malaria care. Finally, planned malaria control activities, including the mass LLIN distribution campaign (implemented between June 2020 and March 2021) and IRS rounds, were successfully implemented despite the pandemic.
In addition, the outbreak was expected to increase fever cases because COVID-19 and malaria both present with fevers. However, community transmission of COVID-19 remained low during the first year following the outbreak; this may explain the fact that this analysis showed no increase in non-malarial visits or proportion of cases with suspected malaria. However, a decrease in the proportion of patients tested for malaria using RDTs was observed. This may be because of global supply chain issues that resulted in local stock-outs of RDTs. In the last year, there was a noted global decline in malaria RDT production and supply, with manufactures shifting in focus from RDTs to COVID-19 point-of-care tests [26,27]. Importantly, the decrease in proportion of RDTs used as a laboratory diagnostic did not affect whether or not patients received a diagnostic test for suspected malaria because microscopy was available as an alternative diagnostic at the MRCs.
There was also a modest decrease in the mean proportion of malaria cases prescribed AL, the standard of care for the treatment of uncomplicated malaria in Uganda, seen in the last six months of observation. Given that a major limitation of interrupted time series analyses is that a change in a time series could be due to another factor that co-occurs with the intervention studied (in this case the arrival of the COVID-19 pandemic in Uganda), it was prudent to investigate whether there were any other changes affecting AL use at health facilities during the COVID-19 time period. Recently, Uganda's Ministry of Health secured funding to procure DP for use as second-line treatment for uncomplicated malaria in the event of treatment failure for patients using AL, and for prescription to all patients with severe malaria following completion of three intravenous artesunate doses [28]. Therefore, from 2020 to-date, DP has been available at some public health facilities for use as a secondline treatment for patients with treatment failure due to AL and all patients diagnosed with severe malaria following completion of artesunate treatment, leading to an increase in the number of patients prescribed DP as described in the results. The decrease in prescriptions of AL observed in the last six months of the study may not be due to a decrease in supply or availability of AL due to the COVID-19 pandemic but instead due to changing prescription practices in response to policy changes regarding DP. There was no direct data on AL stock captured at these health facilities over the timeframe of the analysis, and as such this study was unable to conclusively attribute this decline in AL use to supply chain difficulties caused by the COVID-19 pandemic.
Other studies in similarly resource-restricted settings have shown varied impacts of the COVID-19 pandemic on malaria burden and deaths. A brief report published regarding patients at one health facility in Sierra Leone showed a decrease in the number of malaria cases diagnosed in children younger than five years of age in April 2020 compared to April 2019, but this had recovered by May 2020 [29]. However, this study was at a single health facility, did not account for seasonality, and only observed three months of the COVID-19 time period. A larger study in Zimbabwe, which utilized routine malaria surveillance data collected from all public and private health facilities from January 2017 to June 2020, found an excess of over 30,000 malaria cases from January to June 2020 compared to an average over the same period from 2017, 2018, and 2019 data [30]. This study also found that the number of malaria deaths recorded from January to June 2020 exceeded the annual totals for 2018 and 2019. However, while they describe these increases by comparing numbers from January 2020 to June 2020 to prior data, they do not make any statistical inference regarding their findings and note the limitations of routine surveillance data regarding data quality. A particular strength of the study presented here is the use of a quasi-experimental, interrupted time series study design which accounted for pre-COVID trends, other factors such as seasonality, rainfall, differences in study site, and concurrent LLIN and IRS campaigns, to predict unobserved counterfactual outcomes one year after the onset of COVID-19. In addition, the MRCs included had a full three years of pre-COVID data that had undergone a supervised data cleaning process, resulting in less than 5% missing data.
Several limitations of this study are worth noting. Data from only 17 MRCs was used for this study because highquality data over a sufficient length of pre-intervention time are important to establish a baseline trend in a time series analysis, and most MRCs did not have three years of pre-COVID-19 data. The majority of the MRCs included in the study are in rural areas and, as seen in Fig. 1, are not evenly distributed throughout the country. Therefore, this study's findings may not be generalizable to the entirety of the country, especially to geographic areas not represented in the study and to more urban areas. Though similar outcomes are available through national health management information systems (HMIS) databases, HMIS data have been shown to suffer from increased 'missingness' and quality control problems [31,32] and to provide robust estimates these data were not included. In addition, the impact of the COVID-19 pandemic on malaria mortality was not assessed, as these data were not available. This study was also limited in its ability to assess what drugs patients actually took since data were only available on antimalarial prescription practices and not whether anti-malarial drugs were actually available and administered.
Importantly, COVID-19 cases were still relatively low in Uganda during the period this study was conducted, and it is likely that the impact of the COVID-19 pandemic on malaria indicators may change as community transmission increases, especially with the Delta variant now in the country. This is already being observed as the number of cases between March and July 2021 has doubled compared to what was observed between March and July 2020 (90,656 versus 40,962) [33]. Thus, the conclusions drawn in this study only reflect the effect of the COVID-19 pandemic on malaria indicators in the first year after arrival of COVID-19 in Uganda. The surge in cases caused by the Delta variant resulted in a second government-imposed lockdown in June 2021. Therefore, additional studies on the impact of the epidemic and the effect of lockdown policies will be required.
Finally, though ITSA is known to be a robust method of utilizing observational longitudinal data for estimating causal effects [14,16], the study cannot rule out that estimates may have been subject to unmeasured confounding. For the AL outcome in particular, other health policy changes may be responsible for the decline in prescriptions of AL seen in the second six months of the COVID-19 time period. However, this does not change the central finding, which showed very little effect of the COVID-19 pandemic on malaria burden and case management in the first year after the arrival of COVID-19 in Uganda. Furthermore, if the decline in RDT use is confirmed to be due to supply chain issues and stock-outs, this is critical feedback for the Ministry of Health, so that plans can be made to source additional RDTs and/or to support the availability of microscopy at health facilities.

Conclusion
Utilizing routine surveillance data from 17 rural health facilities across Uganda, this study found the onset of the COVID-19 epidemic did not have major effects on malaria disease burden and indicators of case management, except for a modest decrease in the mean proportion of RDTs used for malaria diagnosis and the mean proportion of malaria cases prescribed AL in the second half of the COVID-19 year. However, continued surveillance will be essential to monitor for any changes in trends in malaria indicators due to the indirect effects of the COVID-19 pandemic so that necessary public health interventions can be instituted early to prevent increases in malaria morbidity and mortality. For example, now that a decrease in the usage of RDTs has been shown, it will be important to further investigate the cause of this decline and implement solutions by sourcing additional RDT supplies or ensuring access to high-quality microscopy at all public health facilities.